The excited state structure of the nitrogen-vacancy center in diamond 
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Optical and microwave double resonance techniques are used to obtain the excited state structure 
of single nitrogen-vacancy centers in diamond. The excited state is an orbital doublet and it is 
shown that it can be split and associated transition strengths varied by external electric fields and 
by strain. A group theoretical model is developed. It gives a good account of the observations and 
contributes to an improved understanding of the electronic structure of the center. The findings are 
important for quantum information processing and other applications of the center. 



PACS numbers: 78.20.Jq, 61.72.Ji, 78.55.-m, 76.30.Mi,42.50.Gy 



Impurity spins in solids are appealing for quantum in- 
formation processing and information storage and the 
nitrogen- vacancy center (NV) in diamond is a defect cen- 
ter with a spin ground state that offers great promise 
HI S IE H ■ It nas been shown that it is possible to opti- 
cally detect the center at single site level Q . Broadband 
optical excitation causes preferentially population of one 
of the ground state levels Q providing a simple way of 
initializing the quantum system. More significantly it has 
been shown that using one laser excitation it is possible 
to read-out the spin state of a single site 0, @ and using 
two-laser excitation to manipulated spins of a single site 
0- Conceivably these capabilities can offer the possibil- 
ity of all optical control of the electron spin. What is 
not clear is whether the processes can be achieved simul- 
taneously with the same NV center and this is vital if 
optical control is to be realized. The readout requires an 
almost perfectly cyclic spin-conserving transition 0, Q 
and the manipulation requires a non-spin conserving or 
A- type transition ^(j- These are contrasting situations 
and the observations may have involved centers in very 
different environments. It is not known how a cyclic or 
A-type transitions is obtained and this is a consequence 
of having a poor understanding of the excited state of the 
NV center. The present work has been directed at im- 
proving this situation. Bichromatic laser and microwave 
optical double resonance measurements are made to ob- 
tain the excited state structure of single NV centers and 
variations with external electric fields are also obtained. 
The observations are compared with a group theoretical 
model and good correspondence is obtained. The model 
provides the necessary insight to the excited state struc- 
ture and enables the immediate questions regarding the 
transition types to be answered and as to how the tran- 
sitions maybe optimized. However, the work has wider 
implication as it leads to an excellent understanding of 
the electronic structure of the NV center and will enable 
more satisfactory development of all NV applications. 

The NV defect in diamond comprises of a nitrogen 



atom at a lattice site next to a carbon vacancy giving 
a center with C3„ symmetry. The center gives an al- 
lowed transition between an orbital A2 ground state and 
an orbital E excited state ^l| . Both the ground and ex- 
cited states are spin triplets (S = 1) 0. The spin levels 
in the 3 A2 ground state are split by 2.88 GHz into a spin 
singlet S 2 and a spin doublet (S x ,S y ). From hole burning 
measurements 0, ^ is known that the 3 E excited 
state is likewise split by several GHz but the details have 
not been determined. The 3 A2 - 3 E line widths are of 
the order of several MHz 0, El and it should be possible 
to establish the excited state structure from an excita- 
tion spectrum. However, excitation spectrum of single 
NV centers have given only one line |7J, la Il3j rather than 
the several anticipated. The reason for the smaller num- 
ber is attributed to differences in the axis of quantization 
in ground and excited state such that with one laser for 
some of the transitions the center is immediately driven 
out of resonance by non-spin-conserving electronic tran- 
sitions. Therefore, the complete 3 A2 - 3 E fine structure 
cannot be obtained with one-laser but can be obtained 
with two and this is the approach used here. 

A laser is modulated to give a beam with a central 
frequency and sidebands at ± 2.88 GHz (Fig 1(a)). The 
laser is tuned to the 3 A2 - 3 E zero-phonon line at 637 nm 
and focused to within a type Ha diamond sample at 4.2K. 
The vibronic emission at 650 - 850 nm is detected in a 
confocal arrangement to obtain the dual-laser excitation 
spectrum of single N-V centers ,13]. Using modulation, 
spectral responses are obtained at several wavelengths 
whereas there are no equivalent responses in the absence 
of the modulation. An illustration of this situation is 
given by the central and lower spectra in Fig 1 (b). The 
laser-sideband separation matches the ground spin split- 
ting and responses are obtained when two of the optical 
frequencies become resonant with transitions from the 
ground state spin levels to a common excited state level 
in a A-type scheme. As the ground state levels are associ- 
ated with distinct spins states, S z and (S^Sy), there will 
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FIG. 1: (a) Energy levels and laser excitation fields. 
(b)Spectral trails of the same NV center when irradiated with 
a laser plus sidebands at ±2.88 GHz (lowest), with single laser 
(middle) and with a laser plus microwaves at 2.88 GHz (high- 
est). 
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FIG. 2: (a) Energy levels and assignment of transitions, (b) 
Spectra taken using single laser plus microwaves at ±2.88 GHz 
(top) and using lasers plus sidebands (second top). The lower 
traces are obtained from simulations including optical and 
microwave fields for a seven level model. The slightly larger 
line widths obtained in the experiment are attributed to laser 
frequency drifts during the experiment. 



only be allowed transitions when the excited state has a 
mixed spin character. The mixing occurs most efficiently 
between adjacent excited states and both of these mixed 
states in general will satisfy the conditions for emission. 
Hence, there will be two features in the sideband (SB) 
excitation spectrum. This is consistent with observation 
and in Fig 1(b) and Fig 2(b) where the two features are 
separated by 0.32 GHz. The features are assigned to the 
1-2 and 3-4 transitions in the energy level scheme in Fig 
2(a). The spectral features are repeated in the SB spec- 
trum at an interval of 2.88 GHz when the laser and the 
reverse sideband are resonant with the transitions. 

In a second approach microwaves (MW) are applied 
resonant with the ground state transition at 2.88 GHz 
to maintain a time-averaged population within each of 
the ground state spin levels. A single laser is then swept 
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FIG. 3: A2 - E excitation spectra as a function of external 
electric field, (a) Two sets of excitation lines associated with 
a single NV center obtained using laser beam modulated at 
2.88 GHz. (b) The variation of the excitation spectra of a 
different NV center as a function of electric field. In both 
figures some of the traces have been manually enhanced. 



in frequency to give a more conventional excitation spec- 
trum. Such a MW spectrum is shown in the upper traces 
in Fig 1(b) and Fig 2(b). The transitions, 1, 2 and 4 to 
the two excited state levels separated by 0.32 GHz are 
observed. An additional transition, 6, is observed and 
interpreted to be to an excited state level at 1.5 GHz 
(Fig 2(a)). 

The measurements in Fig 1 and Fig 2 are interpreted in 
terms of transitions between three spin levels S z , S^, S y in 
both the ground and the excited state each considered to 
be orbital singlets. The interpretation is consistent with 
that given in previous publications |T^. A simulation of 
the spectra using the six spin levels (plus an intermediate 
singlet level) is given in the lower traces in Fig 2(b) and 
gives reasonable agreement with the experiments. There 
is still the need to account for the behavior where one 
excitation line is obtained with a single laser. This line 
is attributed to a different orbital level and, hence, the 
two spectra assigned to the two orbital components of 
a strain-split 3 E excited state. If this interpretation is 
correct all the features should be observable for a sin- 
gle center and this is the situation obtained in Fig 3(a). 
The higher energy features can be obtained with a single 
laser and the lower features only when using a modulated 
laser beam. The features in the upper and lower spectra 
exhibit correlated jumps and have orthogonal polariza- 
tion supporting their assignment to the same defect. The 
modulation used to obtain both spectral types in a single 
scan results in the repetition after 2.88 GHz. 

Electrodes were deposited on the <001> surface of the 
diamond to enable electric fields of up to 3 MV m _1 to 
be applied. This field has components parallel, F z , and 
perpendicular, (F x ,F y ), to the axis of the NV center (x 
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is defined to be in the Cay reflection plane). For the 
latter spectrum (Fig 3(a)) the electric field causes the 
higher energy features to be displaced with a curvature 
to higher energy. This is a characteristic that has been 
reported previously [13j . The lower energy features have 
the reverse behavior with a curvature to lower energy. 
The electric field also causes changes to the fine struc- 
ture and typical variations are shown in Fig 3(b). The 
changes in separations and intensities are consistent with 
two avoided crossings of the spin levels occurring as a 
function of electric field 

The observations can be explained by considering the 
effects on the 3 E state with the Hamiltonian: 

H = Ho + H so + Helec + Hstr 

H gives the energy of the state. H so is the spin orbit 
interaction and H e iec and H s t,- are perturbations due to 
electric fields and strain, respectively. In Cs v symmetry 
the S z spin wave-function transforming as A2 and the 
associated 3 E spin-orbit state has A2<8> E = E symmetry. 
Spins S x , S y transform as rows of an E representation and 
the spin-orbit states have E ® E = Ai, A 2 , E symmetries. 
The four spin-orbit states are affected by spin-orbit: 

H SO = A(L^S 2 + L^S-j. + IjySy) 

First consider the effect of the transverse spin-orbit in- 
teraction. The spin operators S x , S y can in principle lead 
to mixing of spin states but the associated orbital coef- 
ficient A<E| \L xy \ |E> = e is small. Initially this value is 
taken to be zero such that there is no S z - (S x ,S y ) spin 
mixing and the wavefunctions are determined by symme- 
try Axial spin-orbit interaction AL Z S Z will not affect 
the ( 3 E)E state associated with S z spin but will displace 
the states associated with S x ,S y spins (( 3 E)Ai,A2) and 
( 3 E)E by A<E||L X , / ||E> = a and - a, respectively. From 
a six electron (two hole) model a is considered to be pos- 
itive [3. 

The electronic model of the NV center anticipates in- 
termediate x Ai and 1 E singlet levels ^4|. Spin-orbit in- 
teraction with the 1 Ai will cause the ( 3 E)Ai state to be 
displaced to higher energy by b and similarly that with 
the 1 E singlet level can cause one of the ( 3 E)E states 
(lower) to be shifted up in energy by c. This latter shift 
could also arise from spin-spin interaction, a, b and c are 
treated as semi-empirical parameters and are estimated 
from experiments. Values used here are a = 4.4 GHz, b 
= 2 GHz, c = 1 GHz and these values determine the zero 
field energy levels in Fig 4(a). 

The electric field perturbation is given by 

Heiec = F Z D Z + F X D X + F y Dy 

where F$ is the electric field in direction i = z, x, y. Dj 
is the associated orbital operator. The first term F Z D 2 
will simply give rise to linear shift of all levels and not 
considered further. The transverse terms will cause linear 
splittings of the orbitally degenerate E state. The split- 
ting is independent of field direction in xy plane and is 

given by ± dF = ± d ^ 1 {F* + F 2 ) where d= <E||D||E>. 

From previous work d is of order of 6 GHz /MVm -1 [l3|. 
The orbital eigenstates vary with direction of the electric 
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FIG. 4: (a) Energy levels of the excited 3 E excited state as a 
function of the orbital splitting calculated for a field. At 
zero field the levels are determined by the spin-orbit parame- 
ters a, b, and c given in text. The wavefunctions |-E x /( s /)y = 

\E X ) ± \E y ) and \S x / iy /)) = \S X ) ± \S y ). (b) Eigenstates giv- 
ing fraction of S z within each state. In the upper case there 
is only the one state with significant S z component. Special 
features are labelled (i) - (iii). 

field and for a F x field are \E X ) and \E y ), and a F y field 
\E X ) ± \E y ). The variation in the energy of the 3 E spin- 
orbit states for a field F y applied along the y direction is 
calculated and shown in Fig 4(a). The dominant effect is 
the linear splitting of the orbital components. 

In the experiment a quadratic dependence is obtained 
and is due to the presence of strain. This differs from the 
interpretation given in reference [13]. The terms within 
H sir transform as Ai and E irreducible representations, 
the same as those for the electric field. Therefore, the in- 
herent strain has the same effect as a weak electric field 
and the components along z, x and y directions are de- 
noted as f 3 , i x and i y respectively. Using the example of 
a F y field the splitting of the orbital levels are given by 
± d \ff%~+ {f y + F y ) 2 . These quadratic variations are 
consistent with the observations shown in Fig 3(a). The 
minimum separation of the two branches is determined 
by the orthogonal strain field f x = 1 MVm -1 and the 
zero off-set is determined by the parallel strain field i y = 
0.15 MVm" 1 . 

The fine structure is changed by the electric field and in 
the upper branch the spin-orbit separations are reduced 
as the splitting increases (Fig 4(a)). In the lower branch 
the state associated with S z spin changes from being the 
highest in energy at low electric field but crosses to be 
the lowest at high field. However, the non-axial spin- 
orbit interaction previously neglected can cause mixing 
and change crossings to avoided crossings. In the case of 
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a F y field as can occur with the electric field aligned with 
a <100> direction there are two avoided crossings, (i) & 
(ii). Setting A<E||L a; y||E> = e = 0.2 GHz the calculated 
variation shown in Fig 4(a) is in good correspondence 
with the experimental situation of Fig 3(b). 

The eigenfunctions are calculated and the spin projec- 
tions associated with S 2 states are shown in Fig 4(b). 
The lowest trace shows the mixing of the basis states in 
the region of the two avoided crossings ,(i) & (ii). The 
variation of mixing is consistent with the observations in 
the two-laser experiment (Fig 3(b)). The situation with 
the upper branch is totally different. In this case the 
S 2 spin level is robust to spin mixing (Fig 4(b)). For a 
small electric field there is less than 1 in 10 3 component 
of other spins (iii) and the mixing remains small over 
the range of electric fields calculated. Transitions from 
the S z ground level to this upper branch level will give a 
cyclic transition and enable it to be observed in a single 
laser excitation spectrum. For the other upper branch 
states there can be decay via a singlet so that the tran- 
sitions are not cyclic and will not be observed with one 
laser. Also as there is no (S x , S y ) - S z mixing the levels 
will not be observed in the two-laser experiments. It is 
concluded that model gives good correspondence with all 
of the general characteristics observed experimentally. 

The model allows us to comment on other issues. For 
example it is noted that external electric field duplicate 
the effects of naturally occurring strain and the cyclic 
transitions reported for strained diamond will be equiv- 
alent to those for electric field above. The cyclic transi- 
tions are optimal at low but not zero strain (Fig 4(b)) and 
can occur when there is a A transition associated with the 
lower orbital branch. This gives a positive answers to the 
question of whether that cyclic and A transitions can oc- 
cur simultaneously for a single NV center. By varying 
the field or strain one transitions type can be optimized 
but the other transition is degraded. It is also possi- 



ble to anticipate the degree of spin polarization. Near 
an avoid crossing situation there will be almost no spin 
polarization when exciting the lower branch but almost 
total using the upper branch, and perhaps 30% when av- 
eraging both as with vibronic excitation. When using 
vibronic excitation the spin polarization can be enhance 
by work well away from the avoided crossing situation 
as can be achieved with strain or electric field splittings 
greater than 15 GHz. Judging from inhomogeneous line 
widths such splittings may regularly occur and with such 
samples high spin polarizations maybe attainable. Other 
properties for ensembles can be obtained by averaging 
over the strain distribution of NV single centers. 

The present work has provided significant insight into 
the 3 E excited state and when combined with previous 
studies has given us an excellent understanding of the 
electronic structure of the NV center. This is invaluable 
for the development of applications for the NV center but 
in particular for its use in quantum information process- 
ing. An all optical control of single electron and probably 
even single nuclear spins is within reach. This might be of 
particular importance for quantum memory and repeater 
schemes |lq where efficient photon to spin state conver- 
sion is of importance. Fundamental issues and limits can 
be determined and it can be seen when and how the sys- 
tem can be optimized for specific applications. Moreover 
this study demonstrates that for the NV center the rela- 
tive strength of spin flip versus spin allowed optical tran- 
sitions can be tuned by an external control parameter. 
This unusual property might facilitate efficient shelving 
and retrieval of quantum states. 
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